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ABSTRACT 

Theory suggests that a star making a close passage by a supermassive black hole at the center of a 
galaxy can under most circumstances be expected to emit a giant flare of radiation as it is disrupted 
and a portion of the resulting stream of shock-heated stellar debris falls back onto the black hole 
itself. We examine the first results of an ongoing archival survey of galaxy clusters using Chandra and 
XMM- Newton-selected data, and report a likely tidal disruption flare from SDSS J131122. 15-012345. 6 
in Abell 1689. The flare is observed to vary by a factor of > 30 over at least 2 years, to have maximum 
L x (0.3 - 3.0 keV) > 5 x 10 42 erg s -1 and to emit as a blackbody with kT ~ 0.12 keV. From the 
galaxy population as determined by existing studies of the cluster, wc estimate a tidal disruption rate 
of 1.2 x 10~ 4 galaxy -1 year -1 if we assume a contribution to the observable rate from galaxies whose 
range of luminosities corresponds to a central black hole mass (M.) between 10 6 and 10 s Mq. 

Subject headings: X-rays — Galaxies: X-rays — bursts, Galaxies: clusters: individual: A1689, Galax- 
ies: Active, Galaxies: Nuclei — X-rays 



1. INTRODUCTION 

It is now generally accepted that massive black 
hole s (MBHs) reside at t he centers of massive galax- 
ies (|Magorrian et al.lll998|) . These central massive (> 
1O 6 M ) black holes are expected to play an integral role 
in the formation and evolution of galaxies, but many un- 
certainties remain in this picture. 

The low-mass cutoff for the central black hole mass 
spectrum is unknown and the population is inferred from 
the empirical relationship between MBHs and their host 
spheroid' s velocity dispersion and associated luminos- 
ity ("e.g.. iMagorrian et all 119981: iMarconi fc Hunt! 120031: 
lLauer et al.l |2007|) or AGN-ba sed models of growth by 



black hole (jKhokhlov fe Melial fl996l ) . Although much 
theoretical work has been done to describe the expected 
evolution of the debr is and behavior of the resulting flare 
from an analytical (lLoeb fe Ulmerl Il997t lUlmerl 1 1999 . 
e.g.) or numerical approach (lEyans fc Kochanekl 11989 : 
' Khokhlov et~aBl993tlAval et al.l l2000HBogdan ovic et al l 
20041: iStrubbe fc Quataertll2009fl . on the order of a dozen 



accretion ( Marconi et al.|[20(34T ) . While some dwarf galax- 
ies appear to host MBHs and are classified as AGNs, the 
active fraction decreases with decreasing mass and many 
common assumptions regarding the black hole distribu- 
tion with respect to ga laxy mass break down at lower 
M, |G reene et "ail 12008 ?) . Indeed, an unknown fraction 
of dwarf galaxies may host no MBH at all but m erely a 
nuclear stellar cluster (e.g. iFerrarese et al~ll2006| ). 

Outside of observations of the Galactic center 
(|Ghez et alJl2003HGenzel et al.ll2003l ). direct evidence of 
central black holes in the form of accretion events in non- 
AGN (normal) galaxies has been sparse. The most di- 
rect observational evidence of accretion events takes the 
form of flares produce d by a star that passes close to 
the central black hole (jReesi fl988l ) . The passage is so 
close that the star is disrupted by the tidal forces it ex- 
periences. In this case, > 1/2 of its mass is expelled 
into space as the remainder falls back to accrete onto the 
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solid obse rvational examples have been so classified to 
date. See lGezari etal] (l200l for an summary of obser- 
vations and ([Alexan der 2005) for a review of the theory. 

Almost all the observations of tidal flares (or flare 
candidates) have been found in surveys of the field 
(lEscmei et alJ [20071 [20081: iDonlev eFail l200l or in in 
the UV (jGezari et al.l I2009L 120061 ). and while our work 
was already in progress has there been an event reported 
to ta ke place in a cluster of galaxies (jCappelluti et al.l 
2009). 

Our paper is the first in a planned series of studies 
aimed specifically at monitoring clusters of galaxies for 
flares using archival Chandra and XMM-Newton data. 
This paper is based on observations of the first cluster 
we have chosen to study, Abell 1689 (Luminosity dis- 
tance modulus 39.68, z ~ 0.1820) . A1689 was an ideal 
choice for this survey because it exhibits many charac- 
teristics which maximized the chance of observing tidal 
flares given existing rate predictions. 

Abell 1689 is among the richest of the Abell clusters, 
one of only 9 that has a richness of 4 or higher (Abell 
665 is the only cluster of richness 5). At a redshift of 

1 From the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and 
Space Administration. 
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z ~ 0.18, with a virial radius r2oo of 1.13 h^ 1 Mpc and 
hence about 10', the entire core of the cluster fits within a 
single Chandra ACIS chip (|Andersson fc Madeiskll2004j ). 
Its richness has made it popular for examinations of the 
Sunyaev-Zeldovich effect, structure analysis, and lensing 
studies so it is well-observed in X-rays, optical and other 
wavelengths (e.g. Hubble) and hence an ideal starting 
point for our survey. With this extensive coverage (de- 
scribed in detail below) of a rich cluster of galaxies whose 
members are all effectively at the same redshift, we have 
thus been able to estimate the rate at which the galaxies 
experience tidal flare events. 

In this paper we present the discovery of a tran- 
sient event that we demonstrate to be most proba- 
bly due to emission from a tidal flare and whose tem- 
perature is arguably the highest yet observed. For, 
although the derived peak observed flare temperature 
matc hes the highest reported peak blackbody temper- 
ature ICappelluti etHI (pOOll . our result was based on 
data from an instrument that provides about ~ 3 times 
better energy r esolution than was used to derive the tem- 
perature of the ICappelluti et al.l (|2009D flare. 

Furthermore, by monitoring so many galaxies it the 
field of view (FOV) at the same time, it was pos- 
sible for us to set a limit to the frequency of flare 
events in a well defined galaxy popula tion. Depending 
on the hypothesized event rate (e.g., iWang fc Merrittl 
120041: IMerrittl 120091 ) . our results may be used to de- 
limit either the event rate or the very existence of mas- 
sive black holes in dwarf galaxies. Tighter constraints 
on the event rate, as derived from the work presented 
here, are also interesting for predictions of event rates 
of the Light Interferometer Space Antennd^ (hereafter 
LISA) which will be sensitive in the mass range of 
M, < 10 7 (lJennrichll200l ISesana et al.l l200i [20051: 



ISigurdssonl[2003T iKobavashi et al. 2001" 

Throughout this paper, we adopt concordant cos- 
mological parameters of Hq = 70 km -1 sec -1 
Mpc -1 , ^ m .n =0-3 an d £1\.q =0.7, and calcuate dis- 
tances using IWrightl (|2006f ) . All coordinates arc 
J2000. The galactic column density of neutral hydro- 
gen for Abell 1689 is 1.8 2 x 10 20 cm -2 , derived using 
IDickev fc Lockmanl (I199C ) valu es from the colden tool 
in CI AO ( Fruscione et al.lf2006|) unless otherwise stated. 
All X-ray fluxes and luminosities used in this paper have 
been corrected for Galactic absorption using these values. 

2. OBSERVATIONS AND DATA 

2.1. Overview 

In order to identify tidal flare candidates, we begin 
by compiling a catalog of X-ray point sources and iden- 
tify the characteristic properties that distinguish such 
flares from other X-ray point sources such as AGNs, su- 
pernovae (SNe), Gamma Ray Bursts (GRBs) and flar- 
ing M-dwarf stars. We can expect most variable X-ray 
sources to be AGNs. Unlike AG Ns which typically hav e 
a hard power-law spectrum (e.g. iBeckmann et alll2006l ). 
tidal flares will have extremely soft (kT < 0.1 keV) black- 
body spectra. Tidal flares should have much greater vari- 
ability than AGNs near their peak, but that timescale 
may be of short (< 1 yr) duration. A demonstrated 

2 http://lisa.nasa.gov/ 



pre-flare quiescent period as characterized by a lack of 
significant X-ray emission would also distinguish a tidal 
flare from an AGN. Redshiftcd background AGNs, su- 
pernovae, GRBs and M-dwarfs may also have soft X- 
ray emission, but background AGNs may be identified 
by their optical spectra and M-dwarfs will also produce 
characteristic optical emission. 

X-ray-luminous SNe present a subtle contamination 
challenge to an archival survey given the likely impracti- 
cality of optical follow-up on long timescales and possible 
lack of sufficient light curve coverage to describe a decay 
index. But SNe may be eliminated as an explanation 
through a combination of the flare's peak luminosity and 
observed rate of decay. SNe and GRBs typically produce 
X-rays on a timescale of days or months and may be 
be identified optically as well. If the strongest modes of 
soft, luminous X-ray emission arise first from the prompt 
shoc k breakout from the su pernova progenitor's surface 
(e.g. iSoderberg et al.ll2008|) . then from the expansion of 
the ejecta as they interact with the surroun ding medium 
(ISchlegel fc Petrdl2l)0l Hmmler et al.l[2008h . SNe may be 
correspondingly distinguished from tidal flares. 

The m ost X-ray-luminous SNe due to ejecta expan- 
sion (e.g. Ilmmler fc Lewinl[200l ISchlegel fc Petrel 120061: 
llmmler et al.l l2008f) may have long (> years) periods of 
shallow decay. But whether described by thermal or 
power law spectra, they have not been demonstrated to 
have peak luminosities significantly in excess of ~ 10 42 
erg s _1 and typically do not exceed ~ 10 40 erg s _1 . The 
progenitor breakout phase may have a supersoft ther- 
mal X-ray luminosity in excess of 10 43 erg s _1 , as is also 
expected from a tidal flare, but this breakout has also 
shown a fast ri se and exponent i al dec ay on a timescale 
of ~ 100 sec (ISoderberg et all 120081) . Lack of strong 
variability over timescales of typical Chandra or XMM- 
Newton galaxy cluster observations (tens of kiloscconds) 
should thus exclude breakout from the progenitor as an 
explanation for a supersoft X-ray flare of comparable lu- 
minosity. 

The location of the event relative to the host galaxy can 
be used to rule out a tidal flare event if the location is not 
consistent with the core of a host galaxy, such as for a 
supernova outside the galaxy core or a background GRB 
with no host. As such, we have sought to associate flare 
candidates with host galaxies using optical data, and to 
constrain the Rayleigh- Jeans tail of any candidate events 
with UV data where available. 

For our primary dataset, we obtained Chandra ACIS 
archival data for six epochs spanning a period of 8 years, 
from April, 2000 to March, 2007. In addition, we ex- 
amined XMM-Newton data from December, 2001. X- 
ray epochs are indicated in Table [TJ In order to fol- 
low up specific sources, we supplemented our dataset 
with archival HST WFPC2 images with filters F606W 
and F814W. For our candidate event, we obtained a 
March 2009 spectrum of the putative host galaxy with 
the Hobby — Eberly Telescope (HET). 

2.2. Chandra X-ray Observations 
2.2.1. Data Reduction and Candidate Selection 

We reprocessed the ACIS event files using CIA O to re- 
move pixel randomization. According to standard CI AO 
data processing threads, we also applied TGAIN and CTI 
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corrections, and VFAINT mode cleaning where neces- 
sary. We filtered the events to use grades 0, 2, 3, 4 and 
6 and energies between 0.3 and 10 keV. During all obser- 
vations, the target was placed at the ACIS-I aimpoint. 
All observations used ACIS chips 10-4 and S2, except for 
obsid 7701 which used the ACIS-I chips only. ACIS suf- 
fers from charge readout issues which may be remedied 
through use of a destreak software tool, but the prob- 
lem is generally only significant on the S4 chip which 
was inactive in all epochs. We thus did not apply de- 
streak. To account for background cosmic-ray induced 
flares, we filtered individual chips for 3<r deviations us- 
ing the CXC analyzeJtcrv.sl script after excluding bright 
sources found using wavdetect. Minor flaring, of duration 
20 ksec or longer, was present in all epochs; data from the 
corresponding time intervals were excluded from further 
processing. 

We defined several band s in accordance with the 
ChaMP definitions found in Kim et all (|2007t) : B (0.3- 
8 keV), S (0.3-2.5 keV), H (2.5-8 keV), SI (0.3-0.9 keV), 
and S2 (0.9-2.5 keV). In each band, we produced 
weighted exposure maps and compiled a source list with 
wavdetect. wavdetect used a y/2 scale progression of 
1, V2, 2, 2^2, 4, 4^2, 8 on band images binned by a factor 
of two, covering a square of width 4096 pixels in phys- 
ical coordinates that enclosed the whole of the ACIS-I 
array and hence the cluster core. We used a significance 
threshold of 10~ 6 for wavdetect, corresponding to an ap- 
proximate maximum of one false source per chip when 
singly binned, and ~ 1 false source over the whole array 
when binned by two. 

We also created merged, fluxed mosaic images in each 
band with 1" pixels, which we broke into overlapping im- 
ages with dimensions of 1024 x 1024 pixels for the pur- 
poses of running wavdetect. 

Once we had compiled a master source list, we elimi- 
nated conspicuous edge sources and consolidated redun- 
dant multiple detections of single sources across epochs, 
retaining the astromctric positions that demonstrated 
the best accuracy. We also eliminated most sources 
within 0.'5 (~ 95 kpc) of the cluster center, determin- 
ing them to be unreliable based on MARX simulations 
using a f3 model for the Intergalactic Medium. We had 
found these sources to be indistinguishable from Poisson 
fluctuations in the intracluster diffuse emission. We then 
visually selected overlapping sources for consolidation, 
using the 95% encircled energy radii as a guideline and 
giving preference to positions of detections with larger 
SNR, or to Chandra sources in case of additional detec- 
tion by XMM-Newton. 

For initial X-ray s ource classificati o n, we determined 
count rates using the iFeigelson et al.l (|2002f ) polynomial 
approximations for the 95% encircled energy radius of 
each epoch, r 95 and a background annulus between the 
99% encircled energy radius r 9g and 5 x r 99 . We used 
dmextract to determine count rates for each source and 
determined a subset of sources that across epochs demon- 
strated (0.3-2.5 keV) var iability by a f actor of a few, 
> 4 — 5. All errors use the lGehrelsl IpSfih approximation 
for Poisson statistics. 

Within the set of sources that demonstrated signif- 
icant variability, we examined basic spectral proper- 
ties to determine likely soft blackbodies and rule out 



AGNs at low-to-moderate redshifts, which we expect 
to be the most common point source at galactic lati- 
tude b > 60. We define hardness ratios between bands 
= (C2 — c\)/(c2 + Ci) where c„ is the net number 
of counts in C/iaMP-dcfincd band n (minimum c n = 0). 
Given a typical AGN continuum with X-ray photon in- 
dex T < 2.5 (F x oc E~ r ), HR H ,S2 will be sensitive to 
unusually soft spectra except in the presence of a signif- 
icant hard background. Even at Lx(0-3-8 keV) 3> 10 41 
erg s _1 , the exponential cutoff on the Wien side of a 
soft blackbody spectrum will result in few events in the 
H band, if any. But HRs2.si will be an exceptional 
diagnostic for blackbodies in the temperature range an- 
ticipated for tidal disruptions if peak emission is some- 
where in the vicinity of the Chandra bandpass cutoff at 
kT < 0.3 keV. 

We also compared variable sources to the Sloan 
Digital Sky Survey (SPSS) ca talog, Data Release 6 
(jAdelman-McCarthv et aLll2008| ) , eliminating candidates 
that corresponded to a photometric rcdshift consistent 
with background objects (z ^> 0.2), most of which were 
also flagged by SDSS as likely QSOs. We also eliminated 
stars, determined to be bright non-QSO objects desig- 
nated as stars by SDSS. 

Of the 226 distinct sources we examined, almost all 
varied in count rate in one or more X-ray bands across 
epochs, but such variation was usually statistically in- 
significant. For all but four sources, counting statistics 
in the hard (2.5-8 keV) band were insufficient to conclu- 
sively demonstrate variability. Among the highly vari- 
able sources, i.e. those which varied by x4 or more in 
B (0.3-8 keV), only a few demonstrated an unambigu- 
ous early X-ray nondetection (such that the nondetec- 
tion could not be attributed to the choice of the back- 
ground extraction region and local background fluctua- 
tions), followed by single follow-on rise in photon flux, 
as well as later flux evolution consistent with steady de- 
cline across epochs. Of this subset of sources, only one 
(source 141) found at (a, 6) = 13 ,l ll m 22 s , -1°23'45" 
(J2000) both demonstrated an unusually soft peak spec- 
trum and was not eliminated via SDSS as a back- 
ground object, a likely QSO or a foreground star. All 
other highly variable sources were also associated with 
AGNs/QSOs or foreground stars, and also had harder 
spectra (HRh,s > —1.00). Except for sources 78 and 
166, whic h are associated with SDSS stars that have 
NOMAD (jZacharias et al.l 120051) proper motions > 14 
mas yr -1 , and source 153 whose variability is inconsis- 
tent with a flare's rise and fall, all other sources have 
HRr.s > -0.64. 

The X-ray images of the source observed across all 
epochs and its immediate vicinity are displayed in Figure 

m 

For comparison, Table [2] shows a subset containing the 
22 most significantly variable sources, such that all had 
count rates which varied by both at least 2er and by x4 
over at least 2 epochs of Chandra coverage. Due to vari- 
ations in aimpoint, roll angle and choice of active ACIS 
chipset across epochs, many sources fell within the Chan- 
dra FOV for a fraction of the available observations, as 
indicated in Table [2] ■ Rough broadband (B) count rates 
and hardness ratios (HR = HRh,s) were derived using 
dmextract with a constant extraction radius of 10", suf- 
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ficient to enclose 3> 90% enclosed energy even at large 
aimpoint offsets, and a concentric background annulus 
with twice the area of the extraction radius. Note that 
for small 9 this will introduce additional rate uncertainty, 
due to the modest size of the background annulus and 
due to the enclosure of additional background within a 
larger source extraction radius. But it also improves the 
tolerance for positional uncertainty which may be an is- 
sue for fainter sources. 

In Figure [2j we have plotted HRh,s vs. peak count 
rates from Table [2] Compared against the subset of 
significantly variable sources not corresponding to stars, 
source 141 is a clear soft outlier. Most of these non-stellar 
sources cluster at HRh.s ~ —0.50, but source 141 is the 
lowest. 

In very few cases is it possible to place strong limits 
on the pre-peak absence of an X-ray source. Peak emis- 
sion in the first epoch of coverage accounts for 8 of these 
objects. In most cases where the observed peak is later 
than the first epoch of coverage, manual examination of 
the event image provides at least some evidence of pre- 
peak source emission (except for the three objects noted 
in Table [2]) . With those sources eliminated from consid- 
eration, all remaining significantly variable sources could 
not have had steady emission stronger than the la de- 
tection level of 1.9 x 10~ 4 ct s _1 , determined by a single 
10 ks observation. 

Most objects have one corresponding SDSS object 
within the Chandra PSF (a few "). But almost all of 
these objects have been flagged due to traits that are 
problematic to interpretation as a tidal disruption. SDSS 
has flagged 12 optical matches as likely stars due to ob- 
ject morphology and psf modeling, of which 4 are also 
flagged as potential QSOs due to colors atypical of stars. 
If these flags are accurate descriptions of the objects in 
question, they would tend to exclude tidal disruptions in 
normal galaxies as an explanation for the X-ray variabil- 
ity. Exceptions to this restriction can be imagined, such 
as if the SDSS observation were near a flare's X-ray peak 
and captured the blue Raylcigh- Jeans tail in the case 
of an object flagged as a possible QSO. Bright emission 
from the center of a relatively faint galaxy might also 
mimic a stellar radial profile. But without strong evi- 
dence from the X-ray emission, the simpler explanation 
of a star or QSO is preferred. Matches can also be found 
between 4 of these star-li ke sources with SDSS correspon - 
dences and the 1RXH (IROSAT Scientific Team! 12000ft 
and 2RXP (jROSAT Consortiumll200'o¥ catalogs, suggest- 
ing persistent pre-existing sources rather than new flares. 

Excluding all 12 of these flagged objects from Tabled 
leaves 10 sources, of which only 3 significantly variable 
objects are classified as SDSS galaxies. SDSS photomet- 
ric redshifts indicate source 167 is a likely foreground 
object, source 157 is a likely background object, and 141 
is a plausible cluster member. The 7 other unflagged ob- 
jects have no corresponding SDSS objects and therefore 
two possible interpretations. Either they must be back- 
ground objects which we would not examine as flares 
for the purpose of this paper (and for which the most 
likely interpretation is an AGN, given the compositio n 
of the X-ray backgro und, see e.g. iGiacconi et al.l (|2001l ); 
ISzokolv et al.l (|2004f ) and references therein) , or as dwarf 
galaxies below the optical detection threshold, for which 
a tidal flare would remain a possibility of interest. Two 



such sources, 134 and 153, appear to have an optical 
counterpart unregistered by SDSS due to proximity to 
an image flaw. But for the unmatched sources, the X-ray 
lightcurves are either too faint or too sparsely sampled 
to make strong claims as to whether they are the result 
of tidal disruptions. Including those tagged as possible 
QSOs, there are 7 star-like SDSS-matched sources which 
do not have X-ray variability demonstrably inconsistent 
with a tidal flare. In §4 we will consider those as possible 
contributors to the flare rate, as well as the 7 unmatched 
objects, but the available data do not warrant more de- 
tailed examination. 

The exact behavior of tidal fallback emission remains 
uncertain and may be expected to harden according to 

T bb oc M. 1/4 fe.g. lUlmerlll999l ). But HR indicates an in- 
terestingly soft spectrum from the sole likely A1689 clus- 
ter member that is both definitely an SDSS galaxy and 
significantly X-ray variable. According to HR, it is one 
of the two major soft outliers in the significant variability 
population (the other of which, source 153, demonstrates 
erratic light curve behavior clearly inconsistent with the 
simple rise and fall of a tidal flare). 

Having excluded the other X-ray sources from consid- 
eration as tidal flare candidates, we will examine source 
141 in detail in the remainder of §2 and §3. 

2.2.2. Photometry of Source 141 

High State: — The source had been observed in ob- 
sid 5004 on February 28, 2004. It had been found 
by wavdetect in this Chandra epoch only, and there 
in all bands except H. Within rgg, Chandra detected 
a total of 41 counts in the X-ray B band over 19.9 
ksec, with an exceptionally soft HRh,S2 = —1-00 and 
HRs2,si = —0.49. Within r 95 this becomes 37 total 
counts and HRs2,si = —0.41. With an estimated r g5 
background contribution of 5.45 counts, this is a ~ 9.0a 
detection above background fluctuations with a net rate 
of 1.6 ± 0.3 x 10~ 3 s" 1 . Given the H band source emis- 
sion is indistinguishable from background, this becomes 
~ 11.4(7 when restricted to the S band only. 

At 3'7 from the cluster core, contributions from the dif- 
fuse ICM and faint galactic haloes are sufficiently signif- 
icant that care must be taken in choosing a background 
extraction region. Along the radial direction from the 
cluster core the background varies by a factor of > 2 
within 20" of the source for an X-ray B band image of 
all epochs merged. But the source is sufficiently bright 
in the 2004 epoch that choosing a background annulus 
within 5 x rgg does not significantly affect source esti- 
mates. 

Pre-flare: — With no detections of the source in two 
previous Chandra epochs of only 10.5 ksec apiece, we 
coadded results from obsids 540 and 1663 using epoch- 
appropriate extractions to determine more stringent up- 
per bounds for a source that is assumed not to be sig- 
nificantly variable. Respectively, the rgg count rates are 
within 0.5cr or less than those of the background extrac- 
tion annulus. There arc no X-ray B band counts within 
rso for either epoch. Over both pre-flare epochs we find 
a count rate of 9.5 x 10 -5 s -1 within rgs, still less than 
the estimated rgs background rate of 2.5 x 10~ 4 s" 1 . We 
therefore estimate a la upper bound of 1.3 x 10~ 4 s _1 
in excess of the background. 
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Decay: — A much fainter source was also found by 
wavdetect at the same coordinates after the 2004 out- 
burst, but only in the S and SI bands of obsid 7289 be- 
ginning March 9th, 2006. Over 76 ksec there are 51 total 
and 11.33 net X-ray B band counts within r 99 . Using 
smaller extraction radii results in the loss of a small frac- 
tion of the source counts but also reduces contributions 
due to background, which generally increases with in- 
creasing radii due to an increased fraction of the brighter 
core regions of the ICM. Thus, in determining the signifi- 
cance of the March 2006 emission, which is faint with re- 
spect to the background, we use the 90% encircled energy 
radius enclosed by a background annulus out to 2 x rgo. 
From the March 9th, 2006 observations, we find 28 total 
and 11.48 net X-ray B band counts within 76 ks, i.e., 
~ 2.2<t above background. Although there is no detec- 
tion from obsid 6930 on March 6th, 2006, we similarly 
find 18 total and 6.98 net counts within 77 ks, i.e., ~ 1.7<7 
above background. The two March 2006 observations are 
separated by a time interval equivalent to their combined 
duration, i.e., much smaller than the ~ 2 yr elapsed from 
obsid 5004. Because of this short separation, we hence- 
forth treat obsids 6930 and 7289 as a single 153 ks ob- 
servation, which leads to a ~ 3cr detection of the source 
with a count rate of 1.2 ± 0.5 x 10~ 4 s -1 . 

At only 5 ksec, obsid 7701 produced no detection of 
the source and no photons within rgs on March 7, 2007. 
Taking the emission to be background, we estimate the 
la upper bound of 3.7 x 10 -4 s _1 in excess of the back- 
ground. 

Variability: — In order to confirm the degree of variabil- 
ity, we used the CI AO tool glvary, which characterizes 
a Bayesian analysis of Poisson statistics with a variabil- 
ity in dex (VI) based on an logiq(O) where O is the odds 
ratio (|Gregory fc Loredol[T992T lRotil200^ . Examining 
both the X-ray S and B bands across all epochs and 
across expected rise and fall times, we find VI> 8 and 
O > 5.9 in all cases, where VI > 6 is considered definitely 
variable and the maximum VI is 10. The source shows no 
significant variability within its high-state epoch, obsid 
5004. 

2.2.3. Least Squares Spectral Fitting 

We fitted the spectra from the high and post- 
decay states (from February 2004 an d March 2006 , 
respectively ) usin g XSPEC vl2.4.0 (|Arnaudl H995 
iDorman fc Arnaudl 120011 ) with a variety of blackbody 
and power law models that incl ude absorption by line - 
of-sight neutral hydrogen^ and iChurazov et al.1 ()1996l ) 
weighting. We assumed that the source was a member 
of the cluster at fixed z ~ 0.183, although z did not 
change by more than a few percent when left free to vary. 
For spectral fitting we used 99% encircled energy extrac- 
tion radii and account for Chandra ACIS instrumental 
response by calculating position-dependent response ma- 
trices using mkrmf and mkarf. 

We found the high state emission fitted well a black- 
body of temperature of kT^b = 0.12 ± 0.02 keV, with 

3 Due to the relatively low redshift of the object, it does not 
matter if the absorption is assumed to be in our galaxy or the 
external galaxy. For simplicity we assume the absorption all takes 
place within our galaxy for the purposes of spectral fitting. 



X 2 l v = 5.89/5.00. Fixing kTbb at this value, we esti- 
mate an X-ray flux. F x (0.3-10 keV)= 3.7 ± 1.4 x 10' 14 
erg cm' 2 s' 1 , with negligible contribution above 2 keV. 
Fitting tended to d ecrease estimated Nh belo w the galac- 
tic values found bv lDickev fe Lockmanl (|1990l ) to unphys- 
ically low numbers, so we froze Nh at the colden value 
of 1.82 x 10 20 cm -2 . Corrected for galactic extinction, 
F x (0.3-10 kcV)= 4.7 ± 1.7 x 10' 14 erg cm' 2 s" 1 . 

The fit for a power law produced a photon index T = 
4.86 ± 0.48 with a modestly better fit, x 2 l v = 5.36/5.00 
and uncorrected F x (0.3-10 keV)= 5.6 ± 2.5 x 10' 14 
erg cm -2 s' 1 , 99.5% of which is below 2 keV. Allowing 
Nh to vary freely increases it to 1.07 ± 0.41 x 10 21 cm' 2 
with T = 5.8 ± 3.5, x 2 /v = 4.98/5.00 and uncorrected 
F x (0.3-10 kcV)= 5.6 ± 2.5 x 10' 14 erg cm' 2 s _1 . 

For the decayed state, blackbody fitting produces la 
upper limits of Nh = 1.6 x 10 22 cm -2 and kTbb = 0.17 
keV with x 2 l v = 8.35/9.00, consistent with high state 
properties, and a model kTbb = 0.10 ± 0.02, if we fix Nh 
during the fit at the observed galactic value. If we assume 
no change in spectral shape over 2 years, we estimate an 
uncorrected F x (0.3-10 keV)= 4.1±2.8 x 10' 15 erg cm' 2 
s' 1 . Corrected for galactic extinction, Fx(0.3-10 keV)= 
5.1 ± 3.5 x 10' 15 erg cm' 2 s' 1 . 

A power law fit for the post-decay state with galactic 
N H gives r = 5.7 ± 1.1, x 2 lv = 8.93/10.00 and uncor- 
rected Fx(0.3-10 keV) < 1.4 x 10' 14 erg cm" 2 s _1 . T 
cannot be reasonably constrained, if Nh is left free to 
vary. 

An overlaid spectrum of a simultaneous blackbody fit 
of the high and low states is presented in Figure [3] If we 
assume the combined source flux from pre-flare epochs 
represents an approximately constant background, we 
can fit it to high state blackbody temperature, Nh and 
redshift. This fit predicts no pre-flare source Fx (0.3- 
10 keV) with a la upper bound 6.9 x 10' 15 erg cm' 2 
s' 1 . Corrected for galactic extinction, Fx (0.3-10 keV) 
< 8.6 x 10 -15 erg cm -2 s _1 . If we use rgo to reduce 
likely background contamination of this estimate, we find 
an extinction-corrected limit of 1.6 x 10' 15 erg cm' 2 
s . With a much shorter observation time than any 
epoch, we can constrain the post-flare flux from the 
March 2007 epoch only very roughly. Using similar meth- 
ods to those by which we estimate the pre-flare flux, 
even counts within rgs lead to a blackbody extinction- 
corrected Fx (0.3-10 keV)= 8.3 x 10" 15 erg cm" 2 s" 1 
or a la upper bound of 3.5 x 10' 14 erg cm' 2 s' 1 for 4 
counts within rgg in a single energy bin. 

2.3. Supplementary Observations 
2.3.1. XMM-Newton Data 

Pipeline processing for XMM-Newton found no corre- 
sponding X-ray source for the 40 ksec December 2001 
epoch. Our examination of the MOS and PN event 
files and the composite broadband image (0.2-12.0 keV) 
found no evidence of a source within XMM-Newton rgo, 
which at ~ 20" is significantly larger than the Chandra 
positional uncertainty. The mean value of the pipeline 
sensitivity map within rgo is 0.0022 s _1 , equal to the 3cr 
upper limit of the encircled photons. We fitted a PN 
spectrum generated using the espec tool from the XMM- 
Newton SAS software suite, version 7.1.0. For a la black- 
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body with peak kT^, Nh and redshift, we find an uncor- 
rected upper limit (0.3-10 keV)= 8.5x 10" 16 erg cm" 2 
s" 1 or extinction-corrected Fx (0.3-10 keV)= 1.1 x 10" 15 
erg cm -2 s _1 . This source region is also included with 
the Chandra data in Figure [TJ 

We also examined optical and ultraviolet images ob- 
served using the XMM-Newton Optical Monitor with the 
filters UVW2 (1800A-2300A), UVW1 (2300A-3700A), 
U (3000A-3900A), and B (3800A-4900A). The region 
within ~ 2' of the source is dominated by a ghost ring in 
the Optical Monitor, an artifact of increased background 
due to reflection of stray light from a bright source. This 
contamination dominates across all filters and prevents 
extraction of useful optical or ultraviolet data. 

2.3.2. SDSS and CTIO Data 

We examined the Sl oan Digital Sky Survey Catalog 
from Data Release 6 (|Adelman-McCarthv et all 12003 1 
to eliminate flaring red dwarfs as a possible explana- 
tion for variable, extremely soft emission within our field 
and found several such objects corresponding to oth- 
erwise interesting X-ray candidate sources. From our 
softest sources we found one SDSS object designated a 
galaxy by PSF modeling that corresponded to the Chan- 
dra wavdetect position of SDSS J131122. 15-012345.6. 
This galaxy was observed on March 3, 2000 and found 
to have SDSS ugriz magnitudes of u = 25.27 ± 1.14, 
g = 21.45 ± 0.07, r = 20.46 ± 0.04, i = 20.10 ± 0.05 
z = 19.63 ± 0.11. The large u error is likely due to be- 
ing a nondetection in the u-band, which we confirmed by 
examination of ugriz band images. Note that the char- 
acteristic limiting magnitude for SDSS in the it-band is 
22.3 so the uncertainty is likely even higher. 

SDSS Photometric redshifts vary by ~ 0.23 depending 
upon the method used but suggest the galaxy is likely 
a member of A1689 given the much smaller odds of a 
foreground nonmember and typical unreliability of blind 
photometric redshifts at z <C 1. The SDSS template fit- 
ting algorithm finds a redshift of z = 0.097±0.03 whereas 
the catalog's neural network derived redshifts are z = 
0.20 ± 0.09 and z = 0.33 ± 0.07 for the CC2 and ZD1 
methods respectively. We have also run our own pho- 
tometric redshift estimates, finding z = 0.16 ± 0.09 us- 
ing A NNz neural networking software (|Collister &; Lahavl 
l200l . 

We estimate d the redshift with LePhare template fit- 
ting software (jllbert et al.1 12006( 1 and its Hubble Deep 
Field North librar ies. By incorp o rating CTIO J- and 
K-band data from IStanford et all (|2002f >. where Mj = 
19.12 ±0.12 and M K = 17.02 ±0.11, we find an effective 
upper bound of z = 0.20 and a most likely model of an 
elliptical galaxy. The redshift probability distribution is 
bifurcated with maxima at z = 0.04 and 0.20. The upper 
maximum reflects an Sa model, which fits z = 0.19±0.04 
if the elliptical model is eliminated from consideration, 
and a bolometric absolute magnitude Mboi = —20.17 at 
that redshift. Redshift z = 0.04, on the other hand, 
would imply Mboi = —16.29. The fit cannot exclude 
the possibility of a QSO at z ~ 3.44 which assumes Lya 
emission that is undetected in the gap between the g and 
r bands. 

2.3.3. HST-WFPC2 



We retrieved archival images from the Hubble Space 
Telescope taken using the Wide Field Planetary Camera. 
SDSS J131122.15-012345.6 was observed as part of Abell 
1689 by WFPC2 for multiple > 900s orbits on March 
2, 1996 and July 4, 1996 using the F606W and F814W 
filters, although a single orbit from each filter was suffi- 
cient for our purposes. To improve relative astrometry, 
we shifted the HST image to align with Chandra coordi- 
nates. We determined the centroid coordinates of other 
bright X-ray sources within the WFPC2 FOV and visu- 
ally identified the centers of corresponding HST objects 
that also had SDSS catalog entries. SDSS objects cor- 
responded reliably to Chandra objects and were much 
more numerous than Chandra objects, so we used ad- 
ditional SDSS objects to confirm HST correspondence. 
We find the centroid of the Chandra source to be within 
a Chandra pixel (0'.'5) of the host object's center, but 
the primary source of positional error is the statistical 
uncertainty in the centroid of the X-ray source, ~ l'/8. 
The correspondence of the galaxy to the Chandra X-ray 
source can be seen in Figure [4] Later we will assume the 
X-ray source is indeed associated with SDSS J131122.15- 
012345.6 unless explicitly stated otherwise. 

We confirmed the host object galaxy whose mor- 
phology appears plausibly elliptical but exhibits some 
spiral symmetry, possibly containing spiral arms that are 
either broad or too smooth to definitively demarcate at 
the gal axy distance. Using IDL use r-contributed pro- 
grams ([Landsman! Il993t |Barth| |2 001|) for aper ture pho- 
tometry and fitting the flux to a iSersid (|1963|1 i? 1 /™ ra- 
dial profile, we found a Sersic index n ~ 1.2, which sug- 
gests dominan c e by a n exponential disk model ( see, e.g. 
iBlanton et all ()2003[ ): iGraham fc Driverl (|2005l l. The 
flux in excess of an n — 1 profile lies within ~ 0'.'15 
of the center (a 3x3 grid of detector pixels) and suggests 
a maximum bulgc-to-total luminosity ratio B/T= 0.18. 
The radial profile is a poor fit to a deVaucoleurs profile 
where n = 4, even excluding the core excess, suggesting 
the galaxy is not elliptical. We further confirm a spiral 
disk-and-bulge model by examining the ratio of F606W 
to F814W as a function of radius. The ratio is approx- 
imately flat as a function of radius except within the 
region of core excess, where it demonstrates a sudden 
red spike as might be produced by a spiral galaxy with a 
red core and as opposed to a more gradual reddening in 
an elliptical. We conclude, therefore, the galaxy is most 
likely a spiral. 

2.3.4. HET 

We observed SDSS J131122. 15-012345.6 using the 
Hobby-Eberly Telescope (HET) for a total of 5900 sec- 
onds split into 4 sub-exposures. The first two were taken 
on March 24, 2009 and the other two on March 27, 2009, 
well after we would expect to see any emission lines due 
to the tidal flare itself. We used the LRS (low-resolution 
spectrograph) with the Gl grism (300 1/mm) and the 2" 
slit. This gave us a resolution of 14. 8A or 800 km/s at 
6000 A. The spectrum covered the range 4080 - 10800 A, 
with the peak sensitivity between 5500 and 6000 A. The 
signal-to- noise ratio (hereafter S/N) in this spectrum is 
approximately 4 at 5000 A, 7 at 5300 A, approximately 
flat at a value of 10 between 5500 and 7200 A, and drops 
to values of 3-6 between 7500 and 8500 A. At the longest 
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wavelengths, the noise is dominated by imperfect sub- 
traction of the strong OH emission lines in the spectrum 
of the sky background. 

Absolute spectrophotometry with the HET can suffer 
from large errors because a variable fraction of the light 
incident on the primary mirror is focused onto the sec- 
ondary as it tracks the object across the filed of view. 
Therefore, we have adjusted the absolute flux scale by in 
order to match the SDSS photometry. After diving the 
original flux scale by 2.75 we achieve a very good agree- 
ment with the r, i, and z magnitudes. The observed 
spectrum is shown in Figure [5] By comparison with the 
spectrum of the nearby elliptical galaxy NGC 3379, We 
identify features that we interpret to be the Na I D in- 
terstellar double near 7000A and the Ca II H and K 
lines near 4650A at target-redshifted wavelengths, as 
well as Mg I b absorption near 6170A. Thus, we find 
z = 0.195, in good agreement with photometric results 
from LePhare. These features are marked and labeled in 
Figure [S] The signal-to- noise ratio is insufficient to elim- 
inate ambiguity from the redshift estimate, in part due 
to the lack of strong emission lines. This redshift is sup- 
ported by template fitting with the EZ to ol that was used 
in the VIMOS VLT Deep Survey fVVDS. lLeFevre et al.l 
120051 ) and includes spectral line identification capabili- 
ties. The best fits are for elliptical templates which also 
confirm absorption in Ca II and Mg I b, as well as the 
G band, H 7 and H/3 near 5139A, 5186A and 5809A. A 
starburst model also reproduces this redshift but with a 
continuum significantly bluer than observed. None of the 
EZ AGN models converged at the redshift of the object 
(V. Le Brun 2009, private communication). 

In order to evaluate the possibility that emission lines 
form an AGN are hiding in the optical spectrum, we 
have determined upper limits on the flux of some of the 
relevant lines. In particular, we have focused on the ob- 
served wavelength range 5800-6000 A, where we expect 
to find the H/3 and [O III] AA4959,5007 emission lines 
(marked an labeled in Figure [5]). This also happens to 
be the region of the spectrum where the S/N is high- 
est. Although the Ha line is stronger than H/3, it falls 
at an observed wavelength of 7843 A, where the noise is 
very high. The positive spikes discernible in this part of 
the spectrum are a result of poor subtraction of emis- 
sion lines from the night sky. We expect the [O III] 
lines to be unresolved, therefore, we assume a FWHM of 
800 km s" 1 according to the spectral resolution attained 
in the spectrum. Under this assumption we obtain an 
upper limit to the observed flux of the [O III] A5007 line 
of 1 x 10~ 16 erg s _1 cm~ 2 . If a line is stronger than this 
limit, its peak stands out by more than 3tr above the lo- 
cal noise. In the case of the H/3 line, we considered two 
different values of the FWHM, 900 km s _1 (marginally 
resolved) and 1600 km s _1 . We chose these particu- 
lar values because they are typical of the H/3 lines of 
Narrow-Line Seyfert 1 galaxies (hereafter NLSls), which 
are known for their l arge-amplitude X-ray variabilit y 
fe.g. lBoller eTaflll993l . fl997t IGruoe et^ll995all2T)08ah . 
Thus the limits we obtain under these assumptions will 
be relevant to our discussion of whether the observed 
flare could be an example of dramatic variability in a 
NLS1 galaxy. For a FWHM of 900 km s" 1 the upper 
limit to the H/3 is 1 x 10 -16 erg s _1 cm" 2 (just as for the 



[O III] A5007 line), while for a FWHM of 1600 km s" 1 
the upper limit to the H/3 is 2 x 10~ 16 erg s _1 cm~ 2 . For 
line broader than this, the upper limit to the integrated 
flux increases proportionally to the FWHM. 

2.3.5. GALEX 

Ultraviolet observations of SDSS J131122. 15-012345.6 
using the Galaxy Evolution Explore^ {GALEX) exist for 
the dates of 22 April 2004 and 24 April 2007. These 
~ 100s survey observations have limiting magnitudes of 
- 20 - 22 for bands about - 1500A and - 2300A. No 
emission is obviously visible by inspection at the location 
of the galaxy in either band, whether from the galaxy as 
a whole or a point within the galaxy. But given the lim- 
iting magnitude of these short observations, they should 
be insufficient to detect even an unabsorbed 0.12 keV 
blackbody emitting < 4.7 x 10~ 14 erg cm -2 s _1 , as de- 
rived from the Chandra spectrum of 28 February 2004. 
Furthermore, the galaxy lies at ~ 3' of either field's edge 
and is thus flagged for instrumental issues due to rim 
proximity. 

3. DISCUSSION 
3.1. Observational and Theoretical Background 

In order to place our new results in context, we first 
provide brief reviews of both previous observations and 
estimates of tidal flare rates. 

Although we expect tidal disruption flares to be among 
the most luminous observable astrophysical events, with 
the total kinetic energy of ejected debris exceeding that 
of supernovae at 10 51 erg or more, these flares have thus 
far been challenging to observe. While they should oc- 
cur within AGNs and may contribute significantly to 
the faint end of the AG N X-ray luminosity function 
(|Milosavlievic et aLll2006l ). they will be easiest to identify 
in a quiescent galaxy, where they can be distinguished 
from typical variable disk accretion in an AGN. 

Existing theoretical studies predict disruption rates 
of 1 event per 10 4 — 10 5 years per galaxy, (e.g 
iMagorrian fe Tremaind I1999D . a rate that has been 
supported by obs ervational studies using ROSAT 
(jDonlev et al.l 120021 ) . The most optimistic predic- 
tions increase that ra te by an order of magnitude 
(jWang fc MerrittJ[200l . 

Several candidate eyents h a ve been observed (jKomossal 
[200l iKomossa etahl 12001: lHalpern et all I2004D . but 
some of the most convincing evidence of tidal disruptions 
comes from GALEX detections of UV flar es with opti- 
cal and X-ray com ponent as observed by iGezari et al.1 
(poM I2008bl l2009lj. Ongoing studies using the XMM 
Slew Survey (|Esquei et al.ll2007l 120081) also report candi- 
date e vents, but an ext ensive study of the Chandra Deep 
Field (|Luo et al.l I2008D made no detections, consistent 
with maximum rates comparable to 10~ 4 per galaxy per 
year for L > 10 43 erg s _1 . Since only about a dozen 
such candidate events have been identified to date (see 
the above references), the statistical conclusions that can 
be reached thus far are highly tentative. 

The rate by which tidal flares occur should also act 
as an indicator for the distribution of black holes in 
the galaxy population. The effect may be particularly 

4 http://www.galex.caltech.edu 
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pronounced acco r ding to the theoretical calculations of 
IWang fc Merrittl (|200l if a significant fraction of nucle- 
ated dwarf spheroidal galaxies harbor MBHs at their cen- 
ters. Given that dwarf ellipticals are a very numerous 
component of the g alaxy distribution (see, for example 
I Jenkins et al.ll2007f l. if lower mass MBHs flare more of- 
ten than more massive MBHs, they may dominate the 
fla re rate if they c ontribute at all. But more recent work 
by iMerrittl (|2009h suggests lower mass MBHs may pro- 
duce such flares more rarely even if dwarf galaxies do 
possess MBHs. As noted in the introduction, determin- 
ing the population of MBHs in dwarf spheroids through 
such indicators as tidal disruption events will also af- 
fect predicted rates of MBH-MBH mergers and extreme 
mass ratio inspirals (EMRIs) . Although the gravitational 
wave detector LISA will be sensi tive only to disrup- 
tions of stars with compact cores (jAmaro-Seoane et al.1 
120071 ) , as indicators of the dwarf galaxy MBH population 
these predicted rates will have additional implications for 
LISA, which will be sensitive to deeply penetr ating en- 
counters in the mass range of M, < 10 7 M^ ( Jennricbj 
120011 ISesana et alJ[200l 120051 ISigurdssonll2003h . 



3.2. Derived Galaxy Properties 
3.2.1. Cluster Membership 

With a spectroscopic estimate of z — 0.195 for SDSS 
J131122. 15-012345. 6, we take the galaxy to be a mem- 
ber of Abell 1689. A variety of rcdshift estimates 
have been published fo r Abell 168 9 ranging from 0.18 3 
(IStruble fc Rood! [l999h to 0.203 (lEstrada et al.l [2007h . 
The high velocity dispersion of this massive cluster sug- 
gests that membership of the SDSS galaxy is plau- 
sible even for a low estimate of the cluster redshift. 
lHalkola et al.l (|2006t ). for example, estimate a velocity 
dispersion a ~ 1500 km s -1 , placing the galaxy in 
the u pper ~ 1% of clu ster members for the low z esti- 
mate. IDuc et al~l (|2002t ) define an inclusive cluster range 
of z = 0.17 — 0.22 usin g both photometric and spec- 
troscopic redshifts. And iMartini et all (|2007l ) estimate 
a = 2402 km s" 1 for z = 0.1867. For the purpose of 
calculating luminosities then we select a characteristic 
rcdshift z = 0.19 for which the luminosity error will 
be a factor of less than 1.09 whether or not the SDSS 
galaxy's higher measured redshift signifies Hubble expan- 
sion rather than motion relative to the cluster center. In 
our adopted cosmology, z = 0.19 corresponds to a lumi- 
nosity distance of 926 Mpc. 

3.2.2. Black Hole Mass 

Since we do not have a measurement of the stellar 
velocity dispersion of the host galaxy, we have used 
photometry to estimate the mass of its central black 
hole (assu ming it has o ne). We estimated M, from 
My as per ILauer et all (120071) . as well as from Mb as 
per IFerrarese fc FordI (|2005l ). We determine My and 
Mb from SDSS uqriz magnit u des a nd the conversion 
method prescribed by iLuptonl (|2005l ). We specifically 
use the calculations that rely upon the g and r bands 
only given the likely unreliability of u. We have not ap- 
plied K-corrections as the adjustments they introduce at 
z = 0.19 will be minor compared to the errors resulting 
from scatter in the photometric M. — L relation. Using 
the ratio of bulge luminosity to total galaxy luminosity 



(B/T) based on the SDSS data, we find for the bulge 
My = -17.10 ± 0.10 and M g = -15.98 ± 10 

We then re fer to ILauer et all (|2007l ) and 

IFerrarese fc For J (|2005| ) to estimate M, . Given 
My = -17.10, we use the ILauer et all (|2007| ) rela- 
tion that includes their entire dataset, and is hence 
appropriate for My > —19. This relation implies 
lo£{M./M Q ) = 6.23 ± 0.29, and the IFerrarese fctbrdl 
(|2005l ) relation produces log(M./M Q ) = 6.68 ± 0.24. 
Taking the mean from these relations we find 
log(M./M Q ) = 6.46 ± 0.38, an estimate that is 
comfortably compatible with the possibility of tidal 
disruption events, even if a much larger error were 
inferred due to the inherent uncertainty in estimating 
(B/T). The indicated error arises primarily from scatter 
in the relations we reference, but our ability to measure 
(B/T) is limited by WFPC2 pixel scale. For example, 
the host galaxy could have a more compact bulge and 
hence a lower inferred black hole mass. We therefore 
determine the largest uncertainty in the lower bound of 
M, by assuming that the bulge is described by a single 
pixel. In this case, the additional uncertainty due to 
(B/T) is Alog(M./M Q ) cx log{B/T) ~ 0.4. 

3.3. Was the Event a Tidal Flare? 

In the following subsections we argue that alternative 
explanations of the observations are disfavored by the 
data. We then demonstrate that the data are consistent 
with the event being a tidal flare. 

3.3.1. A Galactic Foreground Object? 

At a galactic latitude of 61°, the position of Abell 1689 
already produces a relatively low probability of having a 
line-of-sight foreground object e.g. a distant M-star. We 
now consider the following: (1) A source within the galac- 
tic plane, at a distance of ~ 1 kpc, would have Lx{0.5- 
2.5 keV) < 10 30 erg s _1 , whereas known quiescent low- 
mass X-ray binaries (qLMXBs) typically have Lx{0.5- 
2.5 keV) > 10 31 e rg s -1 and hard spectra, r ~ 1 — 2 
(jHeinke et a,l.ll2003h . (2) M-dwarfs are considerably more 
common and occasionally produce X-ray flares, but even 
a faint galactic M-dwarf should be visible against our 
proposed host galaxy for the event in HST image (take, 
for example, the M-dwarf selec tion criteria of r < 19.0 
by SEGUE, (|Yannv et alJl2009l) ). 

3.3.2. A Highly Variable AGN? 

We consider alternate extragalactic explanations for 
the X-ray flare due to known types of X-ray emitting 
objects. For example, AGNs may show X-ray variabil- 
ity of a factor of a few but typically not bright out- 
bursts of > 30. Characterization of the X-ray continuum 
and the observed power law of T > 5 is also atypical 
of AGNs, which are b etter described by T < 2.5 (e.g. 
iBeckmann et all 120061) . save for NLSls. Although pho- 
tometric fitting is consistent with a QSO at z ~ 3.4, 
luminous AGNs typically have prominent emission lines 
at optical rest wavelengths whereas the host shows no 
strong emission lines in the HET spectrum with limits 
as reported in §2.3.4. A template fit requiring a pho- 
tometric redshift z ~ 3.4 is also contraindicatcd by the 
HET redshift, z = 0.195. Given the optical spectral in- 
formation, we judge therefore that an AGN is unlikely to 
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have caused of the observed flare, with the following two 
exceptions: BL Lac objects and Narrow-Line Seyfert 1 
galaxies (NLSls, known for their steep soft X-ray spec- 
tra and rapid X-ray variability). We consider these two 
possibilities in turn. 

A BL Lac Object. - There exi sts no correspond ing 
radio source in the FIRST catalog (jWhite et allll997D as 
would be expected for a BL Lac. (For a review of BL 
Lac properties, see lUrrv fe Padovani|[l995t) . The FIRST 
catalog reaches a sensitivity 0.75 mJy at 1.4 GHz and 
includes the Chandra fields. Given that most BL Lac 
objects have elliptical galaxies as hosts, the favored spiral 
classification of the host galaxy also implies the flare does 
not arise from a BL Lac object. 

A Narrow-Line Seyfert 1 Galaxy or Similar Object? - 
There are a number of examples of NLS1 galaxies dis- 
playing dramatic X-ray variability (with amplitudes of 
more than a factor of 10). These include WPVS 007, 
which showed a dro p in its soft X-ray f lux by a factor 
of - 400 in 3 years (jGrupe et al.lll995al) . and Mrk 335, 
which showed a drop in its soft X -ray flux by a facto r 
of ~ 30 in approximately 1.4 years (jGrupe et al.ll2008al ). 
X-ray variability by a factor of up to 10 on time scales 
as short as days is a trademark of NLSl s (see, for exam- 
ple, the light curves of NGC 4051 in iMcHardv et~aT1 
l2004f ). An extreme example of this typ e of behav- 
ior is the variability of IRAS 13224-3809 (jBoller et al.l 
fl993l[l997h . There are also examples of dramatic X-ray 
flares in AGNs that are not NLSls, including IC 3599 
(also known as RE J1237+264 or Zwicky 159.034) and 
NGC 5905. Both of these are Seyfert 2 galaxies with 
X-ra y excursions by factors of 120 and 150 respectively 
(see JGrupe et allll995bl; IBrandt et al.lll995l: iBade et all 
fl99l lK omossa fc Badelll999l : IGezari et al.H2003D . 

The characteristic broadband spectral properties of 
SDSS J131122. 15-012345 are generally not unheard of 
among NLSls, but they are atypical in several respects. 
Like NLSls in general, its X-ray emission is very soft. 
The vast majority of NLS1 X-ray spectra, however, ca n 
be fit with power laws such that T ~ 2 — 4 (|Grupdl2004l ). 
Both high and low states of SDSS J131122. 15-012345.6 
are fit to power laws with T ~ 5, which is cqmpara - 
ble to the softest examples from IBrandt et all (|1995l ) 
and moderately steeper than these examples of dramatic 
Seyfert flares, with the dramatic exception of WPVS 007 
(r ~ 10) which h as been attributed to the presence of 
a warm absorber (jGrupe et alll2008b( ). T he correlation 
between spectral slope and Lx observed bv lGrupel (|2004l ) 
is more problematic for classifying SDSS J131122.15- 
012345 as an NLS1. The softest object in that sam- 
ple corresponds to Lx ~ 10 44 erg s _1 , ~ 4x brighter 
than the observed high state power-law Lx (0.2-2.0 keV) 
of SDSS J131122. 15-012345. This high state L x ( 0.2-2.0 
keV) conversely corresponds to T ~ 2.5 in the IGrupd 
(|200l) sample, compared to the observed r ~ 5. The low 
state r is similarly steep, consistent with no time evolu- 
tion, and more difficult to reconcile with this T — L \ 
relationship. Examining the properties of the optical 
continuum, we find that the SDSS g-band absolute mag- 
nitude M„ = —18 .4 is low for the sample examined by 
I Zhou et"all (|2006f ). who find that soft optically selected 
NLSls almost exclusively occupy —19 < My < —25. 
This M g corresponds to a HET spectrum that declines 



steadily with decreasing A for < 5500A (SDSS color 
(u-g)~ 3.8), which would make a similar NLS1 rare 
and " UV -deficient" compa red to a typical flat-spectrum 
NLS1 in (jZhou et al.ll200fih . 

To investigate whether the observed flare in 
SDSS J131122. 15-012345.6 is consistent with the 
variability of an AGN, we ask whether optical emission 
lines characteristic of an AGN could be hidden in 
its optical spectrum below our detection limit. To 
answer this question, we have studied the relation 
between the 0.2-2.0 keV X-ray luminosity and the 
(broad) H/3 and (narrow) [O III] A5007 line lumi- 
no sities among t he 51 NLS1 galaxies in the sample 
of iGrupe et al.l (|2004[ ). In this sample, the ratio 
Fx- (0.2 - 2 keV)/F(H/3) has a median value of 127, with 
90% of objects having F x (0.2 - 2keV)/F(H/3) > 50. 
Similarly the ratio Fxr(0.2 - 2keV)/F([0 III]) has 
a median value of 486. with 90% of the objects 
having Fx (0.2 - 2 keV)/F([0 III]) > 100 Using the 
best-fitting power-law model for the spectrum of 
SDSS J131122. 15-012345.6, we find that in the qui- 
escence F x (0.2-2keV) < 5.1 x 10~ 15 erg cm -2 s" 1 , 
implying that F(H/3) < 1 x 10~ 16 erg cm~ 2 s _1 and 
F([OIII]) < 5 x 10~ 17 ergcm- 2 s _1 (at 90% con- 
fidence), which fall just below our detection limit. 
However, we doubt that this object can be a NLS1 
galaxy because the upper limit on its quiescent X-ray 
luminosity suggests that any AGN should be accreting 
at a rate well below the Edding ton limit, c o ntrary 
to NLSls. In particular, using the lElvis et all (|1994f ) 
spectral energy distribution as a template, the upper 
limit on the quiescent X-ray luminosity, and the black 
hole mass derived above, we infer an upper limit on 
the Eddington ratio in quiescence of -/jboi/^Edd < 0.01. 
In contrast NLS1 galax ies have iboi/^Edd > 0.1 (e.g. 
iMathur fe Grupd [2005). Considering now the cases 
of the Seyfert 2 galaxies, IC 3599 and NGC 5905, we 
notice that in their high X-ray flux states they had 
Fx (0.2 - 2keV)/F([OIII]) = 240 and 10 respectively. If 
these values are representative and if SDSS J131122.15- 
012345.6 is a similar object, we would expect it to have 
F([OIII]) « (1 - 30) x 10" 15 ergem^s- 1 (since it 
has Fx(0.2-2kcV) = 2.7 x 10~ 13 erg cm~ 2 s" 1 for a 
power law in its high state). These line fluxes are an 
order of magnitude above our detection limit, yet we 
have not detected any lines, which leads us to conclude 
that SDSS J131122. 15-012345.6 is unlike IC 3599 and 
NGC 5905. Rathe r it r esembles RX J1242. 6-1119 
( Komossa fc Greinerl [l999h and RX ,11624.9+7554 
( Grupe et alJ 11999( 1 which arc apparently normal 
galaxies that exhibited large, soft X-ray flares with am- 
plitudes of 20 and 235, respectively. Of course, without 
additional observational constraints near the flare peak 
we cannot entirely exclude the possibility that, despite 
having excluded known modes of AGN emission and 
variability we have found an extremely rare instance of 
supersoft variability from a quiescent or obscured AGN 
and not a stellar tidal disruption. Explanations might 
include a novel mode of fluctuation in an otherwise 
low-luminosity (perhaps radiatively inefficient) accretion 
flow or a temporary gap in obscuring material. 

3.3.3. Other Extragalactic Line-of-Sight Objects? 
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As described in §2.1, supernovae have been ob- 
served rarely to emit prompt X-ray flashes that 
reach luminosities comparable to those observed in 
the candi date flare, such as fr o m a progenitor shock 
breakout (ICampana et~aTl 120061: ISoderberg et al.l 120081: 
IGezari et al.l 12008a ). In all cases, the soft thermal com- 
ponent that dominates this prompt flash evolves rapidly 
over timescales of less than the 20 ks peak observation 
of our candidate flare. This evolution appears to occur 
in conjunction with a softening of fcT&fc as the breakout 
expands into larger distances from the progenitor sur- 
face. The A1689 flare demonstrates no such short-term 
variability during the peak epoch, excluding a progen- 
itor breakout as an explanation. The flare is also or- 
ders of magnitude more luminous than even more lu- 
minous exampl es of later, long-term evolution of su- 



pernova eiecta (ITmmler fc Lewinl [20031 iSchlegel fc Petrel 
120061 : llmmler et al.ll2008l ). suggesting that a different ex- 
planation is also necessary here. 

X-ray afterglows from gamma ray b ursts (GRBs) have 
been observed on time scales of years (IGrupe et aLlhoiOfl 
but the decay is sufficiently rapid (declining a factor of 
10 5 within 10 5 s of peak X-ray luminosity, then another 
factor of 10 5 by 10 s s) and the X-ray spectrum sufficiently 
hard (T ~ 1) that we reject them as explanations for the 
event we detected. 

3.3.4. Tidal Flare Explanation 

We now compare our result to previously reported 
events. Early Chandra observations of SDSS J131122.15- 
012345.6 indicate no detectable source for observations 
of at epochs 2000.29 and 2001.98 with L x (0.3-3 keV) 
< 1.6 x 10 41 erg s" 1 at z = 0.19. The XMM-Newton 
observation places a similar limit of 1.1 x 10 41 erg s _1 at 
epoch 2001.98. Then, we detected a prominent ~ IIct 
super-soft flare with an extinction-corrected Lx (0.3-3 
keV) = 4.8 x 10 42 erg s" 1 at t = 2004.16, as well as evi- 
dence (~ 3cr) of emission two years later at t = 2006.18 
with L x (0.3-3 keV) = 5.2xl0 41 ergs- 1 . At* = 2007.18, 
Lx(0.3-3 keV) < 3.6 x 10 42 erg s _1 . Unless otherwise in- 
dicated, further values of Lx for this flare will use a value 
as determined from a 0.12 keV blackbody in the 0.3-3 
keV band, corrected for extinction through the plane of 
our galaxy as per §2.2.3. 

The peak observed X-ray luminosity and temperature 
are as predicted by theory and typical for X-ray obser- 
vation s of tidal flare candidates to date (see IGezari et al.l 
(|2009l) for a review and refer to their Figure 10), and 
the rise in luminosity by a factor of ~ 30 over pre-flarc 
background i s typic al but not remarkable. For example, 
lEsquei et al.1 (|2008f l conside r candidat e s with (0.2-2 
keV) variable by x20 and iLuo et al.l (|2008| ) require a 
count rate variable by > 20, but we see a factor of > 650 
in L x (0.3-2.4 keV) for ICappelluti et al.l <f2009h . 

The observed Lx never approaches th eoretical pea k 
luminosities of > 10 45 for a tidal flare (jUlmerl H999I ) . 
However, the flare may have been observed significantly 
past its peak or could be the result of a "weak" en- 
counter that stripped the outer layers of the stellar at- 
mosphere. Our ability to estimate the flare's bolomct- 
ric luminosity (Ltoi) is also limited by the spectral re- 
sponse of ACIS-I below ~ 0.5 keV. If the spectrum is 
better described as a steep absorbed power law or multi- 



temperature blackbody, Lb i could be an order of mag- 
nitude or more greater than Lx- 

We thus find a tidal disruption event located at the cen- 
ter of SDSS J131122.15-012345.6 to be the most likely 
explanation for the X-ray flare and the following discus- 
sion proceeds based on this assumption. 

3.4. Comparison of Observations to Tidal Flare Models 

With only a few data points we cannot determine the 
shape of the flare light curve observationally. But by ex- 
amining the flare's X-ray light curve and spectral prop- 
erties, we show here that the flare is consistent with the 
expected behavior of a tidal disruption event. 

3.4.1. Light Curve Decay 

Comparison with a Tidal Flare Accretion Model: — At the 
early times, the mass infall rate M for fallback of the stel- 
lar debris from a tidal disruption may exceed the Edding- 
ton limit MeMi even as the luminosit y flare remain s near 
the Eddington luminosity LEdd (e.g. lUlnierl 1 1999ft . The 
central UV/X-ray source may also be obscured at the ear- 
liest times by optically lum inous super-Eddingtqn out - 
flows on a timcscalc of days ([Strubbe fc Quataertll2009( ) . 

At later times, when M < Msddj the luminosity is ex- 
pected to track the rate of return of stellar debris and 
the accretion rate onto the black hole with a decline ac- 
cord ing to \{t — tn ) /(to — tp)1~" where n = 5/3 by the- 
ory (|Phinneylll98ft [Evans fc Kochaneldll98ft lAyal et al.1 
l2000f ). to is the peak M and to is the time of disruption. 

The observations constrain to between the XMM- 
Newton nondctection at 2001.98 and the Chandra high 
state at 2004.16. Even though there are several archival 
observations of Abcll 1689, the flare was clearly detected 
at only only two epochs and was barely detectable by 
2006.18. Given the short duration and the resulting low 
level of statistical significance of the 2007.18 observation, 
we omitted it as a constraint for the purposes of fitting 
the decay of the light curve and simply solve for tu an- 
alytically with two points. For a time t in years where 
t > to > to and to ~ 2003.40, we obtain 



L x (0.3 - 3 keV) = (3.0 ± 2.9) x 10 42 erg s" 1 
~t- (2003.40±0.49) 1_5/3 



1 yr 



(1) 



This relation is shown in Figure [5] with data points, 
error bars and upper limits included. The expected 
rise from to to to is not indicated but is discussed be- 
low. Note the large uncertainty in the normalization is 
strongly correlated with the uncertainty in tjj . 

The variability and spectral shape (see §2.2.2, §2.2.3, 
and Figs. [5] and [3]) of the X-ray data are therefore con- 
sistent with the interpretation that we have identified a 
tidal disruption event, which we will examine in more 
detail as follows. 

Dynamical Interpretation: — By examining plausible val- 
ues for the dynamical timescale of the disruption, we 
place constraints on both the time of disruption and the 
properties of the disrupted star. Theoretically, the dif- 
ference between to and to should be governed by the 
specific energy of the stellar material falling back to the 
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black hole, such as is calculated by iLi et al.l (|2002f ) as 
follows. Let i?„ and M* be the radius and mass of the 
disrupted star, R p is its periastron and Rt is the tidal 
radius of the MBH with mass M.. Then, 

(t - t D ) = 0.11 yr fc" 3 / 2 ( A /' ) V2 

We have increased the timescale h ere according tored- 
shift of the host galaxy, z, as per iGezari et al.l (|2009f ). 
The value of k depends on the extent of the star's spin- 
up by the tidal encounter, with a minimum of = 1 
for a non-rotating star and maximum k = 3 for a star 
that is spun up near the point of break-up. Simulations 
and linear analysis favor values approaching maximum 
value, although the full range is possib le in principle 
([Alexander fc Kumaril200Tl : ILi et al.ll200l . 

If log(M./M Q ) = 6.46 ± 0.38 for the host galaxy, as 
derived from the M. — L relationship, then for a main 
sequence star and k — 3, 

0.02 < (t - t D ) x (R p /R t )- 3 < 0.06 yr. (3) 

This short dynamical timescale implies to -C 2004.16 for 
a main sequence star, suggesting the flare is consistent 
with a tidal disruption event which is already well de- 
cayed from its peak when first observed. 

The analysis above assumes M» ~ M Q and i?* ~ 
Rq, but should also be appropriate for a wide range 
of plausible stellar masses and radii. Note that 
(t -t D ) oc (M./MqJ-^/Rq) 3 / 2 , and while it is 
common to assume a stellar mass-radius relationship 
i?*=R0(M*/MQ) aR wh ere q» = 1 for main-sequence 
stars of M* < 1 M (e.g. lLi et al1l2002D . a R ~ 0.75 is em- 
pirically appropriate for main sequence stars of M* > 1 
i|Hansen et all 120041) . Above M* = 1M , (t - t D ) de- 
pends only weakly on increasing stellar mass. This in- 
cludes the characteristic properties of luminous S-stars 
such as dominate observations of our own galactic center 
(| Alexander! 120051 ) and are therefore of great interest as 
plausible tidal disruption progenitors. 

Suppose an alternative interpretation that does not 
presuppose unobserved early peak emission, such that 
to ~ 2004.16. Such a scenario would be possible if the 
disrupted star had a large radius, such a s a giant, fa- 
cilita ting very long decay timescales (e.g. lLodato et al.l 
2009). The disruption of a giant is plausible given that 
observations of our own galaxy have demonstrated the 
presence of numerous early-type giants, as well as a de- 
pleted popul ation of late-ty pe giants, within ~ 0.1 pc of 
Sag A* (e.g. [Do et aLll200a) . 

In summary, we find that the dynamical time scale of 
the light curve implies a flare significantly decayed from 
its peak luminosity or that the disrupted star was a giant. 
The data do not allow us to distinguish between these 
two possibilities. 

3.4.2. Accreted Mass 

To further examine the properties of the disrupted pro- 
genitor, we estimate the total X-ray energy released by 



integrating the luminosity over the duration of the flare. 
The total mass accreted is thus 

AM = ^| = % [°°L x (t)dt, (4) 

where e ~ 0.1 and fboi is the bolometric correction (here 
~ 1.4). We derived the bolometric correction by tak- 
ing the ratio of luminosities over 0.3-3 keV and over all 
significantly contributing energies as derived through nu- 
merical integration of a 0.12 keV blackbody, for which the 
peak photon energy ~ 0.6 keV. The bolometric correc- 
tion could in principle be an order of magnitude larger, 
however, if instead the X-ray emission arises from a 
steep abs orbed power law or multi-temperature black - 
body (e.g. lGezari et al.ll2009HStrubbe fc Quataertll2009T ) . 

For a lower bound on t, we could conservatively use 
observed luminosity in the high state as the peak of 
flare, which is only ~ 1% of the Eddington Luminos- 
ity L Edd = 1.3 x 1O 45 (M./1O 7 M ) erg s" 1 . This gives 
AM > 1.3 x 10 -3 M Q after the date of peak observed 
luminosity. But if 0.02 < (to — to) < 0.06 yr, as pre- 
dicted from the host galaxy and the M m -a relation, then 
AM ~ 0.01 M Q total, such that M* £ 0.1 M for a 
mass accretion fraction of 0.1. A wide variety of condi- 
tions could produce modest accretion from a more mas- 
sive star. Such scenarios include partial stripping of a 
star's ou ter layers or tidal detonation due to a strong 
counter dCarter fc Luminetl 119821 : iBrassart fc Luminetl 
I2008t lEsquej et all 120081 ) . The range of plausibly ac- 
creted masses is therefore consistent with the mass range 
derived in the previous section. 

3.4.3. Geometry of Emission 

Finally, we use the X-ray spectral properties to address 
the plausibility of an emission region described by a disk 
formed by t he accretion of tidal debris. 

FollowingELelalJ <|2002f >. we can estimate a character- 
istic radius (Rx) of the region emitting X-rays by requir- 
ing the blackbody radiation arise from a projected disk 
and assuming a correction factor f c = 3 to the black- 
body temperature Tm, to account for spectral hardening 
by comptonization. We thus find 

u ( fbol.fc L x \ 1/2 

By this calculation, R x = 1.3 ± 0.4 x 10 12 cm. For 
comparison, the Schwarzschild radius of a black hole 
i?s - 3 x 10 12 (M./10 7 M o ) cm. With only ~ 30 source 
counts, the flare data do not warrant detailed accretion 
disk modeling. But to first approximation we see that 
the area of emission is comparable to the estimated size 
of the host galaxy black hole, Rs ~ 8.6 x 10 11 cm. 

If we assume the X-rays trace the innermost edge of an 
accretion disk, our estimated range of M, from the host 
galaxy also corresponds roughly to a range of innermost 
stable circular orbits (ISCO) for physically permissible 
black hole spins a. such that 0.62 < Rx/Rs ^5 3-6 (com- 
pared to 0.5 < Risco/Rs < 3 for 1 > a > 0, see for 
example iShapiro fc Teukolsky||1986[) . 

3.4.4. Summary of Model-Data Comparisons 
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The spectral properties and light curve of the flare 
are consistent with those predicted for a tidal disrup- 
tion event of a star by a black hole with M, predicted by 
the M,—L relationship. If described by the character- 
istic £ -5 / 3 decay of a tidal flare, the baseline dynamical 
timescale of the galaxy's central black hole implies either 
a disrupted main sequence star that has decayed signif- 
icantly from its peak luminosity by the time of the ob- 
served high state, or a giant star that may instead be near 
to its peak. By integrating the light curve of the flare, 
we find that reasonable assumptions of mass accretion 
rates and efficiency are compatible with the disruption 
of a star where M* ^ 0.1M Q . And a blackbody model 
for the observed spectrum predicts a region of emission 
consistent with an accretion disk, whose inner edge can 
be described by the ISCO over a plausible range of M. 
and black hole spin a. 

4. TIDAL DISRUPTION RATE FROM ABELL 1689 

A comprehensive estimate of the disruption rate per 
galaxy per year is complicated by the nonuniformity of 
the X-ray field coverage, observation spacing, limiting 
flux as a function of position the expected number of 
galaxies within the X-ray fields of view, and the ex- 
pected black hole mass function for different observed 
host galaxy types. A straightforward method of estimat- 
ing the rate is to assume that a given number of events 
N t is characteristic of the number of galaxies observed 
Ng a i and the timespan r over which the observations are 
sensitive to tidal flare emission. Let the rate per galaxy 
per year then simply be 7 = N t /(N ga iT) 

We will first proceed with the conservative assumption 
that the SDSS J131122. 15-012345.6 flare is the only true 
tidal flare. This assumption effectively gives us a "best 
estimate" for the flare rate such that Nt = 1. We can 
estimate the limits of r to which our observations are 
sensitive from the limiting flux of the observations for a 
given black hole mass. Assuming z = 0.19, we can use 
the theoretical decay of flare accretion as described by 
iLi et all (|2002|) to solve for 



it - t D ) 



~L5 



-3/5 



M. 



10 7 M 



3/5 



yr. 



(0) 



where -F-13 is the limiting flux of the observation in terms 
of 10~ 13 erg cm" 2 s" 1 . 

Due to the diffuse X-ray emission of the cluster, the 
limiting flux is dependent on the distance from the clus- 
ter center, where the number of galaxies per unit solid 
angle is greatest. We must thus consider the impact of 
the ICM in our estimation of N ga i and r. To calculate 
our sensitivity, we therefore exclusively use values for the 
S band (0.3-2.5 keV). This selection minimizes the hard 
contribution to the background and we expect few or no 
flare photons at higher energies, whereas the ICM con- 
tribution at higher energies is likely to be significant. 
To determine a comprehensive limiting flux, we select a 
minimum of Lx (0.3-2.5 keV)= 10 42 erg s" 1 for a tidal 
flare high state to differ entiate from Ult ra-Luminous X- 
ray sources (ULXs; e.g.. lHeil et "aTll2009D and other sim- 
ilarly luminous sources, and assume kT^b ~ 0.12 keV 
as observed. This selection criterion also sets a strong 
lower limit to M, for any flares in our sample, given 



L E dd k, 10 42 erg s" 1 for M. ^ 10 4 M Q . The exposure at 
epoch 2006.18 is sufficiently deep that it is sensitive over 
the whole field to disruption events beginning as early 
as the 2001.02 epoch under these selection criteria from 
flares similar to the SDSS J131 122. 15-012345.6 flare such 
that ~ 3cr from M, > 10 6 M Q , except in the central ~ 1' 
of the cluster. There, X-ray emission is dominated by 
the ICM. 

Observations at earlier Chandra epochs (2000.29, 

2001.02) arc sensitive across the field to Lx(0.3 2.5 

keV) ~ 5 x 10 42 erg s" 1 at 3a. Outside of the central 
1' radius, early Chandra epochs are thus sensitive to dis- 
ruptions from M. ~ 10 6 M Q over time scales of ~ 1.5 yr. 
The single XMM-Newton epoch at 2001.98 has a broad- 
band (0.1 — 12.0 keV) 3(7 sensitivity < 10 42 erg s" 1 across 
the field. 

We can estimate N ga i by modeling the number of 
galaxies observed projected in ACIS-I acco rding to pre- 
viousl y determined cluster characteristics. ILemze et al.l 
(|2009D found A1689 to have a limiting radius of 2.1 h" 1 
Mpc assuming z=0.183 and a maximum radial velocity 
amplitude of ~ |4000| km s" 1 at 300/1" 1 kpc. We use the 
surface number density as a function of radius that they 
determined, namely £„ = £o/[l + ( r / r c) 2 ] p + C where 
S ~ 1300 galaxies h~ 2 Mpc" 2 , r c ~ 395 h- 1 kpc, p ~ 1 
and C is the number of background galaxies. We esti- 
mate that a given Chandra field of view with all four 
ACIS-I chips active plus ACIS-S6 contains N gai ~ 2100 
cluster members to a depth of Iab = 26.5, eight magni- 
tudes fainter than the characteristic bend at luminosity 
L* in the cluster luminosity function (e.g. iPaolillo et al.l 
l2001f ). and only ~ 0.1 N ga i will be found in the central 1' 
of bright ICM emission. Since fewer than ~ 10% of the 
galaxies covered by our observations are close enough to 
the cluster core to be effectively obscured by the ICM, 
the bright diffuse background from this region will there- 
fore have only minor effects on the overall estimated rate 
and hence can be ignored. 

If, from the M,-L relation, we assume our candi- 
date flare represents a characteristic M, ~ 3 x 10 6 M Q , 
then the ~ 7 years of X-ray observations for A1689 im- 
ply we are sensitive to tidal flares of similar luminosity 
and duration for at least r ~ 8 years. Combined with 
N ga i ~ 2100, this yields 7 - 6 x 10" 5 galaxy" 1 yr" 1 
for all galaxies in the sample. Without taking effects of 
the galaxy population compositio n into account , this i s 
m excess of - 10" 5 as f ound bv iDonlev et al.l IpOOl . 
within the bounds set bv ILuo et aL (120081) . and abac- 
tor of ~ 4 less than estimate d bv lEsauei et al.l (120081) . 
lEsauei et all (|2008l ) used the iFerrarese fc Ford! ([2005D 
M.-ct relationship and the iWang fc Merrittl (120041) rate 
predictions such that 7 ~ 7.0 (M./lO 6 M )-°- 28 x 10" 4 
galaxy" 1 yr" 1 . For an assumed M, ~ 2.88 x 10 6 M Q , 
the same formulation predicts 7 = 5.2 x 10~ 4 galaxy" 1 
yr" 1 , a factor of 9 greater than the rate we calculated 
above for A1689. 

However, since we expect that lower values of M. will 
have markedly shorter durations based on strength and 
rate of decay, suppose that for M, ~ 3x 10 6 M Q , the 
observable rate will likely be dominated by black holes 
of M, within an order of magnitude. We can estimate 
the contrib uting fracti o n of t he cluster population by in- 
tegrating a iSchechter] (|1976l ) luminosity function. We 
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use K s = 17.02 as above (|Stanford et al.l 120021) and a 
T^-band function as determined by IKing en!T(l2002T) 
such t hat a = —1.01 and to* = 14.6. arcom fc Hunt! 
(|2003h determine an M. — Lk relationship such that 
logipM. oc —2.2K and there is little variation in M m — L 
relationships across several near-infrared bands. We in- 
tegrate the Schechter function over K s ± 2.2, and then 
compare that value to the full galaxy population which 
is integrated down to m± + 8 (an alogously to t he Ia b 
limiting sensitivity determined by iLemze et al.1 (120091) ). 
or M+ + 8 = —16.8 with K-corrections from IKing et al.l 
(|2002j ). We thus find ~ 49% of the galaxies to be likely 
contributors to the production of flares detectable by 
our observations. If, therefore, N ga i sa 0.49 x 2100 « 
1000, the disruption rate we infer from a single event is 
1.2 x 10~ 4 . This r ate is within a factor of two of the 
lEsauei et~a l. (2008) rate that makes no corrections for 
nondetections of flares in heavily absorbed edge-on spi- 
rals. 

But suppose that tidal flares account for some signifi- 
cant fraction of the other sources we detected that both 
displayed significant X-ray variability and could not be 
unambiguously eliminated from consideration. The high- 
est upper bound for the flare rate would include all 14 
objects as potential tidal flares. This bound is too high 
as we can exclude all of the SDSS-matched objects as 
likely AGNs or stars based upon t he optical pr o pertie s 
previously described. According to llvezic et al.l (|2002D . 
SDSS morphological classification is robust to 95% as 
faint as r ~ 21.5, suggesting the matched stellar ob- 
jects (all of comparable magnitude or brighter) are un- 
likely to be misclassified ga laxies. And a tidal optica l 
flare of the type predicted bv lStrubbe fc Quataertl (|2009f ) 
(yL v ~ 10 41 erg s _1 in g) is too faint to be detected by 
SDSS at z ~ 0.19, so that such a flare could not ex- 
plain starlike unresolved optical emission from a cluster 
member. 

Removing these 7 SDSS-matched objects from our 
sample, we are left with 7 remaining unmatched variable 
sources that lack optical identifications and are outside 
the FOV of all HST observations of A1689. These objects 
could be background AGN that are optically fainter than 
the SDSS ugriz limiting magnitudes. We will now calcu- 
late upper limits to the flare rate as determined from two 
alternate cases that do not ignore these unmatched X- 
ray sources. In the first case, we assume that all of these 
variable sources are cluster members. In the second case 
we attempt to estimate what fraction of these unmatched 
sources are line-of-sight objects rather than cluster mem- 
bers. But in either of these cases, any host galaxy 
within t he cluster wou ld have to be a dwarf galaxy. 
With the iLuptonl (120051) magnitude conversions and the 
lLauer et al.l ([2007D M.-cr relation, the SDSS g and r lim- 
iting magnitudes suggest dwarf host galaxies are unlikely 
to have M, > 10 5 5 M Q . A flare without an obvious host 
could also arise fr om an intracluster MBH ej ected during 
an MBH merger (jKomossa fc Merrittl[2008f ). 

In the first case, under the unlikely assumption that 
all optically unidentified sources were indeed tidal dis- 
ruption flares from faint (dwarf) cluster galaxies, the 
rate for all galaxies would be 7 ~ 8.2 x 10~ 4 . But if 
the typical duration of a flare becomes shorter with de- 
creasing M. , the effective value of N ga i r will also be less 



than if we had assumed u niform flare dur ation. Sup- 
pose we therefore assume a iSchechten (|1976l ) profile and 
r cx M 3 / 5 as before, with a power law spectrum T = 2 
(using a harder spectral profile like a typical AGN and 
consistent with the uniform hardness ratio distribution of 
these objects). Rounding our limit of M, < 10 5 5 M up 
to 10 6 M© for dwarf galaxies and given the Eddington- 
limited lower bound of M. ^ 10 4 M Q , the integrated 
N ga iT for 10 6 Mq > M. > 10 4 M Q in our datasct is only 
~ 1400 galaxy-years. Thus, in this case 7 ~ 5.0 x 10~ 3 
for the subset of galaxies where M, < 10 6 M Q . This 
specific rate is several times greater than even the most 
opti mistic values (e.g. 7 = 1 -2 x 10~ 3 for M. = 10 5 M Q 
from lWang fc Merrittl (|2004f )). 

In the second case, we ca lculate a mor e plaus ible dis- 
ruption rate by using the ILemze et al.l (|2009t) galaxy 
distribution function to determine the likelihood of any 
given galaxy being a cluster member rather than a field 
galaxy. This probability decreases at greater radii from 
the cluster center. Starting at r > 800 kpc from the 
cluster center, even the most central of these 7 optically 
unidentified sources would be no more than ~ 40% likely 
to be a cluster member. If we consider only the respec- 
tive projected distances from the cluster center for a set 
of unknown galaxies, there is a 77% chance that at least 
one of these seven galaxies is a cluster member. The most 
likely or "expectation" value is n = 1.3 cluster members. 
Here n = Eo<!<nP( r! ) where n = 7 is the number of 
galaxies and p(r,) is the probability of galaxy i being a 
cluster member at radius r, from the cluster center. The 
expected number of field galaxies is therefore n — h = 5.3. 
This estimate neglects any bias introduced in a dataset 
selected for X-ray emission and variability. 

To determine an upper bound for Nt in this calcu- 
lation, we consider h ~ 1.3 optically unidentified X- 
ray-variable cluster members, plus a tidal flare in SDSS 
J131122. 15-0123 45. 6 and corresp onding statistical error 
of 2.4 based on iGehreld (fl986l ). Therefore Nt < 4.7 
and 7 < 4.8 x 10 -4 . These assumptions also indicate 
7 < 2.6 x 10~ 3 for 10 6 Mq > M. > 10 4 M w , which is con - 
sistcnt with the predictions of IWang fc Merrittl (|2004l ). 
Thus, our observations do not exclude the possibility that 
most dwarf galaxies contain MBHs in this mass range. 

5. SUMMARY AND CONCLUSIONS 

We report an X-ray flare that occurred within a pro- 
jected distance of 1.6 kpc from the center of the spiral 
galaxy SDSS J131122. 15-012345.6. The projected posi- 
tional uncertainty is 4.7 kpc. The peak observed lumi- 
nosity (Lx (0.3-3 keV) ~ 5 x 10 42 erg s _1 , assuming it 
originated from a member of A1689), strong variability 
(x30), supersoft spectrum (kT^ = 0.12 keV), slow evo- 
lution ( x 9 decay over ~ 2 years) , and absence of strong 
optical emission lines from the presumed host galaxy ~ 5 
years after the outburst are consistent with other exist- 
ing examples of X-ray flares that have been classified as 
arising from tidal disruption events. Although we cannot 
completely rule out another origin for the X-ray flare, we 
conclude that the most likely explanation is a flare due to 
the tidal disruption of a star by an MBH. To the best of 
our knowledge, this is the first reported Chandra-selected 
discovery of such a flare. 

Based on observations of a cluster with a well-measured 
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galaxy population and over a well-defined period of time, 
we have therefore calculated an event rate of ~ 10~ 4 
tidal disruption events galaxy" 1 yr _1 under the assump- 
tion that all inactive, non-dwarf galaxies have central 
black holes of M. = 10 6 - 10 8 M and will produce 
tidal flares. The calculated event rate is consistent with 
previously reported theoretical and observational work . 
Our results coupled with those of lCappelluti et af] (|2009f ) 
present a strong case for cluster surveys to detect X-ray 
flares from tidal disruption events. Our work is consis- 
tent with the hypothesis that all normal galaxies (such as 
SDSS J131122. 15-012345.6) host massive black holes at 
their center, but our results do not exclude other models 
of galaxy formation in which MBHs are rarer or (alter- 
nati vely) more massive th an the M, — a relation implies 
(e.g. iVolonteri et "all 120051 ). 

Determining the rate of tidal flares from lower mass 
black holes, if they occur at all, has the potential to 
constrain existing models of dwarf galaxy formation and 
evolution. A more accurate disruption rate for normal 
galaxies in clusters compared to rates derived from the 
field could also be related to galaxy evolution in rich clus- 
ters of galaxies. In particular, it could have implications 
for how galaxy mergers affect black hole mass distribu- 
tion and stellar populations in galaxy clusters. 

Flares from low mass (M. < 10 6 M r] black hole s may 
(|Wang fc Merrittl 120041) or may not (jMerrittl 120091) con- 
tribute disproportionately to the overall flaring rate. But 
if low mass black holes produce tidal flares of shorter ob- 
servable duration, existing cluster X-ray observations of 
a typical ~ 2 year cadence could miss tidal flares even for 
an actual rate that is perhaps a factor of 10 higher. Fur- 
ther dedicated X-ray observations of clusters using Chan- 
dra, XMM-Newton or future soft-sensitive (0.1 — 3.0 keV) 
missions wi th comparable or b etter fi elds of view s uch 
as WFX T (IMurrav et all l2008h . IXO (|Parmarll2009T) or 
eRosita ([Predehl et al.ll2007D are therefore necessary, to 
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CXO: 2000.29 XMM/PN: CXO: 2004.16 CXO: 2006.18 CXO: 2007.18 
CXO: 2001.02 2001.98 

Fig. 1. — Images of the event at the times of observation epochs, starting from pre-fiare. Instruments used are indicated: CXO for 
Chandra X-ray Observatory ACIS, and XMM/PN for the XMM-Newton PN. Numbers correspond to the date of observation. Circles 
describe a 95% encircled energy radius for Chandra ACIS or 50 in the case of XMM-Newton PN. Color scales are uniform and scaled 
to provide a good contrast at epoch 2006.18. Note that this saturates the color scale at epoch 2004.16, the observation nearest to peak. 
Chandra images are binned from (0.3-10 keV) event files by a factor of two, XMM-Newton unfiltered by energy and binned by a factor of 
32. Pre-fiare Chandra epochs have been merged with extraction circles from both epochs indicated (as per differing off-axis angles). 
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Fig. 2. — Sources from Table. [2]with hardness ratios plotted against corresponding peak count rates. The filled circle indicates our tidal 
flare candidate, source 141. Crosses mark those X-ray sources which correspond to SDSS objects designated as stars. Diamonds indicate 
sources all other sources. 
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Fig. 3. — Top: fits to the Chandra spectra using a kT = 0.12 keV blackbody model at the cluster redshift at epochs 2004.16 (solid line: 
model, squares: data) and 2006.18 (dot-dashed line: model, triangles: data). l<r errors are indicated. Each bin signifies 5 photon counts. 
Bottom: residuals between model and data, as above. 
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Fig. 4.— Pre-flare HST WFPC2 image of SDSS J131122.15-012345.6 taken using the F606W filter. The larger overlaid pixels are singly 
binned X-ray events from the flare at is peak at epoch 2004.16. The circle indicates centroid error at r ~1'.'8 and the cross is the middle of 
the centroid. 



A Tidal Disruption Flare in Abell 1689 21 



T 




Ca II 

i i i i I i i i i I i i i i I i i i i I i i 

4000 5000 6000 7000 8000 

Observed Wavelength (A) 

Fig. 5. — Observed HET spectrum of the flaring galaxy with a spectrum of the nearby elliptical galaxy NGC 3379 shifted to z = 0.195 
overlaid for comparison. There is a residual from the imperfect subtraction of the sharp [O I] telluric line at 5577A, marked with a crossed 
circle. In addition, there are also large residuals from poor sky subtraction at wavelengths longer than 7500 A. Amongst the galaxy's 
absorption features, at 7000A we see the Na ID interstellar line, at 4650A are Ca II H and K, and at 6200A is Mg lb. The spectrum near 
Ha is dominated by noise from sky subtraction residuals. Note the lack of strong emission lines that would indicate a typical AGN. 
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Fig. 6. — Light curve for the flare with error bars. Arrows indicate upper bounds. Lx (0.3-3) corrected for galactic extinction are 
indicated by X . The dashed line is a 

£-5/3 i ight 

curve which assumes to at the expected value for a solar- type star, indicated by the 
solid vertical line. Dotted lines indicate the error range for t jj . The rise in luminosity between tr> and to is n °t indicated, nor is the peak 
luminosity at to shown. 



A Tidal Disruption Flare in Abell 1689 



TABLE 1 
X-ray Observations. 







Observation 


Duration 


Aim Point C 


oordinatcs 




Instrument 


Obsid 


Date (UT) 


(ks) 


a(2000) 


5(2000) 


Roll Angle 


ACIS-I 


540 


2000 Apr. 15 


10.5 


13:11:13.88 


-01:20:14.75 


214.302 


ACIS-I 


1663 


2001 Jan. 7 


10.5 


13:11:34.90 


-01:17:10.88 


67.828 


EPIC 


0093030101 


2001 Dec. 24 


39.8 


13:11:33.69 


-01:20:29.70 


115.755 


ACIS-I 


5004 


2004 Feb. 28 


19.9 


13:11:33.30 


-01:17:20.91 


75.154 


ACIS-I 


6930 


2006 Mar. 6 


76.9 


13:11:32.87 


-01:17:21.16 


77.129 


ACIS-I 


7289 


2006 Mar. 9 


75.7 


13:11:32.87 


-01:17:21.17 


77.129 


ACIS-I 


7701 


2007 Mar. 7 


5.1 


13:11:37.04 


-01:19:13.76 


77.454 



TABLE 2 
Variable Sources 



Source 


RA 

{ ) 


Dec 
( J 


SDSS Optical ID 


1 AB 


Chandra Obsid Coverage 
540 1663 5004 6930 7289 7701 


Peak li.atc 
ct s 


Peak HR P 
, H-S 
1 H+S 


'eak Plags 
(') 


14 


197.9353 


-1.1661 


none 






X 


X 


p 


X 




3.18 ± 1.09 


-0.38 


12.8 


1,2 




19 1 .9825 


-1.1709 


J 131155. 5 1-01101 / .( 


18.45 




p 


X 


X 


X 




10.7 ± 4.7 


-0.39 


8.68 


4 


27 


197.8700 


-1.1488 


J131128.78-010856.5 


21.47 




X 


X 


X 


1' 




1.94 ± 1.00 


-0.38 


8.49 


1,3 


31 


197.7475 


-1.1341 


J131059.20-010805.5 


20.77 




X 


X 


X 


p 




5.61 ± 1.47 


-0.52 


12.5 


1,2,4 


34 


197.7035 


-1.1420 


J131049.00-010832.4 


20.15 




X 


X 


X 


p 




5.27 ± 1.46 


-0.42 


14.1 


1,2,3 


13 


197.8888 


-1.2866 


J131133.28-011712.2 


19.71 


P 


X 


X 


X 


X 


X 


34.9 ±7.9 


-0.69 


5.74 


2,3 


75 


197.8232 


-1.2121 


J131117.54-011244.5 


19.92 


X 




X 


1' 


X 


X 


0.26 ± 1.07 


1.00 


6.00 


4 


77 


198.0102 


-1.3321 


none 






X 


p 


X 


X 


X 


6.30 ±2.74 


-0.04 


7.73 


1 


78 


197.9978 


-1.2848 


J131159.55-011705.6 


14.95 




p 


X 


X 


X 


X 


11.7 ± 5.1 


-0.84 


6.15 


2,4 


86 


197.9058 


-1.4036 


none 




X 


X 


X 


X 


p 


X 


4.40 ± 1.43 


-0.52 


6.95 


1 


123 


197.7654 


-1.4008 


none 




P 










X 


15.0 ±5.6 


-0.61 


4.58 


2,6 


131 


197.8080 


-1.4281 


none 




P 










X 


16.9 ± 5.9 


-0.43 


5.44 




134 


197.8259 


-1.2416 


J131118.20-011429.8 


18.36 


X 


X 


X 


X 


X 


p 


277 ± 26 


-0.50 


6.68 


2,3,7,8 


141 


197.8424 


-1.3959 


J131122. 15-012345.6 


i 19.60 


X 


X 


p 


X 


X 


X 


15.6 ±4.1 


-1.00 


6.98 


5 


152 


197.9071 


-1.3304 


J131137.68-011950.0 


20.42 


p 


X 


X 


X 


X 


X 


78.9 ± 12 


-0.64 


5.97 


2,4,8 


153 


197.9144 


-1.3956 


J131139.42-012344.1 


21.29 


p 


X 


X 


X 


X 


X 


5.81 ±5.46 


-1.00 


7.29 


2,4,7 


157 


197.9352 


-1.3191 


J131144.39-011909.3 


20.49 




X 


X 


p 


X 


X 


20.1 ± 20.0 


-0.52 


3.40 


2 


160 


197.9597 


-1.4248 


none 






X 


p 


X 


X 


X 


13.3 ±3.5 


-0.28 


9.19 


2,6 


162 


197.9695 


-1.2442 


none 






X 


X 


p 


X 


X 


9.16 ± 1.52 


-0.35 


5.64 


1,2,5 


163 


197.9758 


-1.3797 


J131154.23-012246.9 


19.99 


X 


p 


X 


X 


X 


X 


23.3 ±6.0 


-0.60 


7.39 


1,2,4,8 


166 


197.9864 


-1.2198 


J131156.76-011311.8 


13.26 




1' 


X 


X 


X 


X 


82.1 ±9.9 


-0.92 


6.76 


2,4,8 


167 


197.9928 


-1.1745 


J131158.32-011026.8 


16.65 




X 


p 


X 


X 




6.54 ± 2.72 


0.62 


9.29 


1,5 



Significantly variable sources and their SDSS optical counterparts, as described in §2.2.1. Bold indicates the flare described in detail in the text. SDSS positions 
may differ slig htly from X-ray source positio ns, but are within X-ray positional uncertainty. Iab magnitudes are calculated from SDSS ugriz magnitudes using 
ILuptonl (|2005T ). Compare to m* = 18.5 from ILemze et all ()2009l ) for Iab- X marks epochs when the object would be located within the ACIS FOV, regardless of 
whether it is detected during the epoch. P indicates the epoch of peak count rate, for which the rate, HR and 6 are presented. H and S bands are as described in the text. 

Flags: (1) High pre-peak uncertainty (2) XMM detected source (3) SDSS flagged as possible QSO (4) SDSS flagged likely star. A small fraction may also be QSOs 
or other bright, unresolved o bjects. (5) Possible pre-peak n ondetectio n (6) SDSS object visible but not catalogued (7) Variability pattern inconsistent with flare (8) 
Pre-existing match in 1RXH (jRQSAT Scientific Teamll2000l ') or 2RXP (|RQSAT Consortium|[2000l ') 



